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The production of prompt photons at high energies provides a direct probe of the dynamics of 
the strong interactions. In particular, one expect that it could be used to constrain the behavior of 
the nuclear gluon distribution in pA and AA collisions. In this letter we investigate the influence 
of nuclear effects in the production of prompt photons and estimate the transverse momentum 



dependence of the nuclear ratios RpA ~ 



dcr{pA) ij^ dtrjpp) 
dyd-^pj' ' dyd'^pj' 



and Raa = 



dcr(AA) I ^2 da{pp) 



dyd-'pT 



/A 



dyd^PT 



at RHIC 



and LHC energies. We demonstrate that the study of these observables can be useful to determine 
the magnitude of the shadowing and antishadowing effects in the nuclear gluon distribution. 
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In the last years the study of heavy ion collisions have 
provided strong evidence for the formation of a quark- 
gluon plasma (QGP) |l|, 0], which is one of the main 
predictions of the Quantum Chromodynamics (QCD). 
Currently, distinct models associated to different assump- 
tions describe reasonably the experimental data, with the 
main uncertainty present in these analysis directly con- 
nected with the poor knowledge of the initial conditions 
of the heavy ion collisions. Theoretically, the early evo- 
lution of these nuclear collisions is governed by the dom- 
inant role of gluons, due to their large interaction proba- 
bility and the large gluonic component in the initial nu- 
clear wave functions. Such extreme condition is expected 
to significantly influence QGP signals and should mod- 
ify the hard probes produced at early times of the heavy 
ion collision. Consequently, a systematic measurement of 
the nuclear gluon distribution is of fundamental interest 
in understanding the parton structure of nuclei and to 
determine the initial conditions of the QGP. Other im- 
portant motivation for the determination of the nuclear 
gluon distributions is that the high density effects ex- 
pected to occur in the high energy limit of QCD should 
be manifest in the modification of the gluon dynamics. 

One of the nuclear effects which is expected to modify 
the behavior of the nuclear gluon distribution is the nu- 
clear shadowing (For a recent review see Ref. Q). This 
effect has been observed in the nuclear structure func- 
tions by different experimental collaborations [J, |5| in 
the study of the deep inelastic lepton scattering (DIS) 
off nuclei. The modifications on F^{x,Q'^) depend on 
the parton momentum fraction x. While for momen- 
tum fractions x < 0.1 (shadowing region) and 0.3 < 
X < 0.7 (EMC region), a depletion is observed in the 
nuclear structure functions, in the intermediate region 
(0.1 < a; < 0.3) it is verified an enhancement known 
as antishadowing. These experimental results strongly 
constrain the behavior of the nuclear quark distribu- 



tions. However, the nuclear gluon distribution is usu- 
ally inferred using the momentum sum rule as constraint 
and/or studying the log Q^ slope of the ratio ¥2"^ /F2 
6]. Due to the scarce experimental data in the small- a; 
region and/or for observables strongly dependent on the 
nuclear gluon distribution, the current status is that its 
behavior is completely undefined. It is demonstrated by 
the analysis of the Fig. [H where we present the results 
for the ratio Ra = xgA/A.xgN predicted by the EKS [3], 
DS 8], HKN d, [i3 and EPS [lli parameterizations at 
Q2 = 2.5 GeV2 and A = 208. These four groups re- 
alize a global analysis of the nuclear experimental data 
using the DGLAP evolution equations [12[. From Fig. [U 
we have that currently the magnitude of the shadowing 
and the presence or not of the antishadowing are open 
questions. It has motivated several authors to propose 
the study of different observables in distinct processes to 
constrain the nuclear gluon distribution (See Ref. [3]). 
For instance, in Ref. [l3| it was proposed to study the 
heavy quark and vector meson production in ultrape- 
ripheral heavy ion collisions at LHC, while in Ref. [1J| 
the study of the longitudinal and charm structure func- 
tions in eA collisions at RHIC. Another possibilities are 
the heavy quark and quarkonium production in central 
proton-nucleus and nucleus-nucleus collisions (See e.g. 
Refs JH [H, [H [S [la, M\- Furthermore, in Refs. 
[21I [23 ] the authors have studied the nuclear effects in 
prompt photon production. This process provides a di- 
rect probe of the dynamics of strong interactions and, 
consequently, of the nuclear gluon distribution, since the 
dominant mechanism of production at high energies is the 
Compton scattering (7 + g ^ 7 + q. More recently, the 
possibility to observe the nuclear gluon shadowing using 
high-pT prompt photon production at RHIC and at LHC 
was discussed in Ref. |2J] . It was argued that the produc- 
tion of isolated photons turns out a promising channel for 
the extraction of Rg . In this letter we extend these pre- 
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FIG. 1: Ratio Rg = xqa/A.xqm predicted by the DS H, EKS 
01, HKN [l^l and EPS [l3| parameterizations at Q^ = 2.5 
GeV'^ and A = 208. The vertical lines indicate the minimum 
X values which are probed in photon production at central 
rapidities in proton-nucleus collisions at RHIC and LHC. 



vious studies considering the more recent nuclear parton 
parameterizations and estimating the transverse momen- 
tum dependence of the ratios R^a — /% I A ^ ''^ and 
Raa = S^M^gl^ at RHIC and LHc'ene'rgies, 



where the differential cross section is given in Eq. ([!]) 
below. Our goal is to determine if these observables can 
be used to constrain the different effects expected to be 
present in the nuclear gluon distribution and to discrimi- 
nate among the distinct models present in the literature. 
Prompt photon production occurs through two types 
of processes: the direct piece, where the photon is emit- 
ted via a point like coupling to a quark, and the frag- 
mentation piece, in which the photon originates from the 
fragnrentation of a final state parton. As the second com- 
ponent can be almost completely reduced by isolation cri- 
terion used in the experimental data analysis, we focus 
our study only in the direct component, which provide a 
clean probe of the hard scattering dynamics. In this case, 
the prompt photon production cross section is given by 

M 



da, 



dydp^ 
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2x\ — XTe-y 



dt 



(Q ,xi,a;2) ,(1) 



where xt = ^pr/y/s, y and pT are the rapidity and trans- 
verse momentum of the produced photon, fi(x,Q^) are 
the parton densities, xi and X2 are the momentum frac- 
tions of the partons involved in the hard process. In this 
case we have that xo = ■^1^^'= " and a;',"™ — ^ — 



2a;i 



2-XTe-y ■ 



Q^ is the hard scale and ^ are the partonic cross sec- 
tions, which are perturbatively calculable [2^. The con- 
tributing LO subprocesses are qg -^ qj (Compton), 



qq — > 57 (annihilation), followed by the subdominant di- 
agrams qq -^ 77 (pure EM), gg — > 77 and gg — > 177. The 
correspondent LO matrix elements and partonic cross 
sections can be found in Refs. [23, |25^- In what follows 
we estimate the differential cross sections for central ra- 
pidities. It implies that the cross section at small values 
of the photon transverse momentum in pPb (dAu) colli- 
sions at LHC (RHIC) is determined by the behavior of 
the nuclear gluon distribution at Xg > 10"'' (10~^) (See 
Fig. [T|). It is important to emphasize that smaller val- 
ues of X contribute at photon production in the forward 
rapidity region which can be measured, e.g., with the 
PHENIX and STAR detectors at RHIC and by the CMS 
experiment at LHC. 

The main input in the calculations of the prompt pho- 
ton cross section are the nuclear parton distribution func- 
tions (nPDF). In the last years several groups has pro- 
posed parameterizations for the nPDF, which are based 
on different assumptions and techniques to perform a 
global fit of different sets of data (in particular, EPS in- 
cludes RHIC data) using the DGLAP evolution equations 
[3,[1,IB,[i3,[11|- In Fig- Uwe present a comparison among 
the distinct parameterizations for the nuclear modifica- 
tion factor for the gluon, Rg{x,Q^), at A — 208 and 
Q^ — 2.5 GeV^. As we can see, these parameterizations 
predict very distinct magnitudes for the nuclear effects. 
For larger values of x, the EKS and the EPS show anti- 
shadowing, while this effect is absent for the HKN and 
EPS parameterizations in the x < 10""' domain. The 
more surprising feature is however the anrount of shadow- 
ing in the different parameterizations. While the shadow- 
ing is moderate for DS and HKN parameterizations and 
somewhat bigger for EKS one, the EPS prediction has a 
much stronger suppression compared with the other pa- 
rameterizations. For smaller x around x ~ 10~^, while 
DS and HKN parameterizations have about 20% suppres- 
sion and EKS one have about 40% suppression, for the 
EPS parameterization this effect goes to almost 80% sup- 
pression in the nuclear gluon compared with the A scaled 
gluon content in the proton! For bigger values of x the 
behavior is distinct for all parameterizations. As x grows, 
the DS parameterization predicts that Rg grows continu- 
ously to 1 , that means that the shadowing dies out when 
X — > 10^^. The same happens for the HKN one in this 
limit, but this growth starts only at a; > 10^^, with Rg 



being flat for 10"^ < a; < 10 



At 



10^ 



we have 



that behaviors predicted by the EKS and EPS parame- 
terizations are similar, with Rg exceeding 1.2. The main 
distinction between these parameterizations is that in the 
EPS parameterization one has a much steeper growth, 
from a much stronger suppression at smaller x to the 
antishadowing behavior for larger values of x. 

In what follows we calculate the prompt photon pro- 
duction cross section pA and AA collisions, considering 
the nuclear parton distributions discussed above, and es- 
timate the ratio between these predictions and the un- 
modified pp production cross section. The so called nu- 
clear modification ratios, defined in the Eqs. (O and ([3]) 
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FIG. 2: Transverse momentum dependence of the ratio RpA 
in prompt photon production at RHIC (y^ — 200 GeV) and 
LHC {^/s — 8.8 Tel^), for different choices of nuclear parton 
distributions: DS ^, EKS ^, HKN O and EPS [H. 



FIG. 3: Transverse momentum dependence of the ratio Raa 
in prompt photon production at RHIC (y^ = 200 GeV) and 
LHC (y^ = 5.5 TeV), for different choices of nuclear parton 
distributions: DS ^, EKS ^, HKN ^ and EPS [11 • 



below, could be measured at RHIC and LHC. We ana- 
lyze the transverse momentum dependence of the nuclear 
modification ratios. One have that if a proton-nucleus 
or nucleus-nucleus collision is nothing more than a su- 
perposition of nucleon-nucleon collisions, the ratios RpA 
and Raa should be unity. As the EKS and EPS sets are 
only evolved to leading order (LO), the calculation of the 
cross sections is also done at LO accuracy. To be consis- 
tent, we pick only the LO version of the other nPDFs we 
are considering. Besides, since we are only calculating 
ratios between cross sections, common uncertainties on 
the normalization of the pA, AA and pp cross sections, 
due to higher order contributions, are expected to cancel 
out in the ratios. 

In Fig. [2] we present our estimates for the transverse 
momentum dependence of the ratio RpA defined by 



_ dajpA) dajpp) 



dydPpT 



dyd^PT 



(2) 



in pA collisions at RHIC (^/i = 200 GeV) and LHC 
(y/s = 8.8 Tel^). The factorization scale Q is assumed to 
be equal to the photon transverse momentum pT- The 
results show distinct behaviors of RpA for the distinct 
nPDF's and for the different colliders. The difference be- 
tween the predictions for RpA at RHIC and LHC comes 
from the kinematical Xg range probed in the two cases. 
While for RHIC the values of Xg are ever larger than 



lo- 



in the LHC case the minimum value is 10 , in- 



creasing with Pt- Consequently, the effects in the nuclear 
gluon distribution which contribute for the prompt pho- 
ton production are different in the two colliders. It is 
verified in Fig. [21 where we observe that at RHIC, the 
EKS and EPS predictions implies RpA larger than one for 
Pt > 5 GeV, being substantially larger for EPS, which 
is directly associated to the magnitude of the antishad- 
owing present in xqa- In contrast, for the DS and HKN 
predictions we have smaller values for the ratio due to 



the absence of this nuclear effect, with the ratio being 
larger for the DS prediction than the HKN one. At large 
Pt, which implies larger values of Xg, the EKS and EPS 
predictions decreases due to the EMC effect. From these 
results we can conclude that the study of the ratio RpA 
at RHIC can be useful to determine the presence or not 
of the antishadowing and constrain its magnitude. On 
the other hand, at LHC energies all parameterizations 
predict a ratio smaller than one, with the suppression di- 
rectly associated to the magnitude of the shadowing ef- 
fect. For instance, we have that the ratio is substantially 
suppressed in the EPS (EKS) case, being smaller than 0.9 
(0.95) at pT < 20 GeV, while in the DS (HKN) case it is 
almost fiat and equal to 0.95 (0.9) in the full pT range. 
Morever, differently from the DS and HKN predictions, 
the EKS and EPS one predict a strong transverse mo- 
mentum dependence. Consequently, the determination 
of the magnitude and pt dependence of the ratio at LHC 
are useful to determine the properties of the shadowing 
in the nuclear gluon distribution. 

The production of photons can also be studied in the 
collision of heavy nuclei. The main interest arises from 
the fact than, once they are produced, photons leave the 
hot and dense fireball virtually without further interac- 
tions (See e.g. Ref. [2^). In particular, prompt photons 
constitute a physics background in the search for thermal 
photons, which are one of the QGP signatures. There- 
fore, a precise determination of the nuclear effects present 
in the initial state of the collision and the consequent sup- 
pression of the cross section is fundamental for the study 
of the photon production in ultrarelativistic heavy ion 
collisions at RHIC and LHC. In what follows we present 
our estimates for the transverse momentum dependence 
of the ratio Raa defined by 
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da {A A) 
dyd^PT 
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(3) 



in AA collisions at RHIC (^/i = 200 GeV) and LHC 
(y^ = 5.5TeV). In this case the nuclear effects are am- 
plified by the presence of two nuclei in the initial state 
of the collision. Similarly to the pA case, we can see 
in Fig. [3] that the ratio Raa is strongly modified by 
the shadowing and antishadowing effects. In particular, 
the ratio is predicted by the EPS parameterization to 
be equal to 1.15 at px — 10 GeV and RHIC energies, 
while the DS parameterization predicts Raa ~ 1.0 in 
the same pr range. Therefore, also in A A collisions at 
RHIC, the study of prompt photon production can be 
used to constrain the magnitude of the antishadowing ef- 
fect. On the other hand, at LHC energies we have that 
the ratio is strongly suppressed by the shadowing effect, 
with Rg being predicted by the EPS parameterization to 
be approximately 0.5 (0.7) at pr = 5.0 (10) GeV . Our 
results indicate that the production of prompt photons 
with pt < 20 GeV can be used to determine the shad- 
owing effect. 

A comment is in order. In this letter we restrict our- 
selves to the descriptions which use the DGLAP evolu- 
tion equations II4I to describe the behavior of the nuclear 
parton distributions. However, at small values of x, it 
is expected the presence of nonlinear corrections in the 
QCD evolution (For a recent review see [23] )■ In this new 
regime, the nuclear gluon distribution should be deter- 
mined by the solution of a nonlinear evolution equation. 



as for example the GLR-MQ evolution equation [28] . As 
shown in Ref. ,29j , where the prompt photon production 
in pp collisions was studied considering the solution of 
GLR-MQ evolution equation for the proton gluon distri- 
bution obtained in [30], new effects are expected when 
nonlinear PDFs are used in the calculations. A similar 
expectation is valid for the nuclear case. 

In summary, in this letter we have investigated the 
prompt photon production in pA and AA collisions at 
RHIC and LHC, considering the coUinear factorization 
and some of the parameterizations for the nuclear par- 
ton distributions available in the literature. Our results 
demonstrate that the nuclear ratios RpA and Raa are 
powerful observables to discriminate among the different 
parton distributions in the nuclear medium. In particu- 
lar, the predicted shadowing by the EPS parameteriza- 
tion is considerably larger than in the previous nuclear 
PDF's. This strong shadowing effect could be tested at 
the forthcoming LHC. 
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